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We present experimental data on the limiting behavior of an interreality system comprising a 
virtual horizontally driven pendulum coupled to its real- world counterpart, where the interaction 
time scale is much shorter than the time scale of the dynamical system. We present experimental 
evidence that if the physical parameters of the simplified virtual system match those of the real 
system within a certain tolerance, there is a transition from an uncorrelated dual reality state to 
a mixed reality state of the system in which the motion of the two pendula is highly correlated. 
The region in parameter space for stable solutions has an Arnold tongue structure for both the 
experimental data and for a numerical simulation. As virtual systems better approximate real ones, 
even weak coupling in other interreality systems may produce sudden changes to mixed reality 
states. 



Although increasingly sophisticated real-time com- 
puter simulations of the real world are created every day, 
to date there has been little to no research done on the 
physics of the pairing of a virtual system and its real- 
world counterpart, often referred to as an "interreality" 
system Virtual systems are often created to model 
real systems as accurately as possible, with great pains 
taken to eke out additional realism. Aside from virtual 
worlds designed for entertainment, examples of accurate 
virtual models of the real world abound in high-energy 
physics accelerator work Furthermore, a computer 
simulation can feature unidirectional coupling to the real 
world, as in the case of black box trading in finance [3] 
or the dynamic clamp in neuroscience [4|. The unidi- 
rectional coupling can also be from the real world to 
the simulation, as in the case of data-driven modeling 
whereby live measurement data is incorporated into an 
executing application. This has been used for creating 
accurate, real-time models of systems ranging from com- 
plex vortex flows Q to human cancer cells 0. The next 
step is to examine an interreality system in which there 
is bidirectional coupling. 

In this Brief Report we present experimental evidence 
for a transition from a dual reality state to a mixed reality 
state in an interreality system featuring bidirectional, in- 
stantaneous coupling. The experimental phase diagram 
is in the form of an Arnold tongue There is good 
agreement between the experimental data and simula- 
tion. We demonstrate that even a simple model taking 
into account a single degree of freedom and using only lin- 
ear damping is sufficient to give rise to this mixed reality 
state. It is often difficult and prohibitively expensive to 
create a model that predicts all observable parameters of 
a system in the real world to maximum precision. Even 
a computer model of the familiar physical pendulum in 
air must take into account linear and quadratic damping. 
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FIG. 1: Experimental apparatus. Fig. [TIa) shows a side view 
detail of the pendulum housing: pendulum (A) is attached to 
slotted disk of angular encoder (B) in housing (C). Fig. [Ub) 
shows a top view of the apparatus: pendulum (A) in housing 
(C) is attached to lever arm (D) with pivot points (E), driven 
by actuator (F). The pendulum has a length of I — 15.37 
cm and the adjustable distance 5.1 cm < a; < 25 cm from 
the pendulum to the pivot point controls the strength of the 
dimensionless coupling constant 0.17 < < 0.67. As indi- 
cated, 6 is the angular displacement of the real pendulum as 
measured from the vertical. 



the buoyancy of the pendulum in the air, the added mass 
due to the pendulum dragging air with it as it moves, and 
a half dozen other effects in order to reproduce the mea- 
sured period to a precision of 10~^ seconds 0. All but 
the most sophisticated computer simulations will use ap- 
proximations when it comes to modeling a physical sys- 
tem. 

As an example of a virtual system coupled to a real one, 
we choose a horizontally driven physical pendulum as the 
real system. The horizontally and vertically driven phys- 
ical pendulum has been described in the literature 
The known equations of motion accurately model the dy- 
namics of the system; to create a physically accurate vir- 
tual pendulum we need only to make a real-time simula- 
tion based on these equations. The experimental physical 
pendulum is a lightweight, very thin wooden rod with a 
length I = 15.37 cm, a diameter of 2.3 mm, and a mass of 
0.4 g. The pendulum is connected directly to the roller 
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of a digital angular encoder. The data from the angu- 
lar encoder is sent to a computer; this provides a simple 
and robust measure of the instantaneous angular position 
of the pendulum. Forcing for the physical pendulum is 
provided by an amplifier and a PASCO actuator (model 
SF-9324) with a maximum displacement of 0.3 cm. The 
pendulum is attached to a lever arm that allows it to be 
driven with a greater amplitude Xdrive , up to a maximum 
of 5.75 cm. The computer controls the voltage to the ac- 
tuator via an analog DAQ. Figs.[T{a) and[IJb) show the 
experimental setup. 

We calculate the equation of motion for a horizontally 
driven physical pendulum in the usual way. See Phelps 
and Hunter for a detailed derivation @. We define 9 
as the angle through which the real pendulum oscillates 
as measured from the vertical, and we define (j) to be 
the corresponding angle for the virtual pendulum. Then 
the equation of motion for the horizontally driven vir- 
tual pendulum with an arbitrary time-dependent driving 
function / is 

-f 2/30 4- (tl'Wo)^ sin0 + a„/„ cos(/i = 0, (1) 

where f3 is the damping coefficient, loq = 2TTUJr is the nat- 
ural angular frequency of the real pendulum, u> = iVv/uJr 
is the dimensionless ratio of the natural frequencies of the 
two pendula. In terms of the pendulum moment of iner- 
tia / and mass m, we define a = I~^mlxdrive x^rive/l 
to be the dimensionless coupling constant that sets the 
scale of the coupling term. Weak coupling corresponds 
to o <C 1. To determine the coefficient of linear damp- 
ing, we fit the decay of the uncoupled real pendulum to a 
simple exponential. For our real pendulum this measured 
value of /3 is 0.45 ± 0.2, while the measured natural fre- 
quency of the real pendulum Ur is 1.57±0.01s~"'^. These 
are the values used for f3 and cOr in each calculation. In 
this work, an overbar on a variable denotes a normalized, 
dimensionless quantity. Also, the subscripts v and r re- 
fer to variables associated with the virtual pendulum and 
the real pendulum, respectively. A standard fifth-order 
Runge-Kutta routine is used to integrate this equation. 
We let / scale with the angular displacement of the real 
pendulum: fv{d) = 0- At the n-th time step of duration 
At, the measured positions of the real pendulum at the 
current and two previous time steps determine the value 
of/: 

- J ^ ^ 

Using these values, the integrator returns </)„, which is 
used to determine the driving amplitude fr^n = 0n for 
the real pendulum. A voltage proportional to /r,n is sent 
to the actuator driving the real pendulum. The distance 
from the pivot point along the lever arm determines 
[see Fig. [TJb)]. Since the measurement, integration, and 
feedback are easily completed by the computer in a time 
Stcomputer with St computer < < I/wq, the program 
then waits to integrate again until the internal timer 



reaches nAt after the initial starting time for the first 
integration. Since it is impractical to attempt to release 
the real pendulum from precisely the same starting po- 
sition each time the program is run, instead the virtual 
pendulum is started with zero initial velocity but with a 
nonzero initial position. The real pendulum is started at 
rest, with 6n=o = 0. Two typical sets of position versus 
time data are plotted in Fig. O 
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FIG. 2: The top plot shows position versus time for the ex- 
perimental data point co = 0.98, d = ^Ja^aZ = 0.19 (Region 
I, mixed reality state). This plot is an example of the system 
exhibiting stable, phase-locked periodic motion. The bottom 
plot has Lo = 0.72 and d = 0.22 (Region II, dual reality state). 
This is an example of the system ending in the stable equilib- 
rium position (p = 6 = 0. For both plots, the solid and dashed 
lines correspond to 9 (the position of the real pendulum) and 
(j) (the position of the virtual pendulum), respectively. 

With At = 35 ms, the bidirectional feedback is ef- 
fectively instantaneous; the effect is that of a real-time 
virtual pendulum simulation that immediately responds 
to any motion of the real pendulum and vice versa. This 
works because At is much shorter than the characteristic 
time scales of the dynamical system. The natural fre- 
quency of the real pendulum is approximately 1.57s~^. 
With no feedback, the motion of either pendulum ceases 
after less than 10 s due to friction in the real pendulum 
and damping in the equation of motion of the virtual 
pendulum. We allow the system to run for 45 s, long 
enough for any transient dynamics to vanish, and then 
we evaluate the final dynamics of the system. We find 
that there are two equilibrium states of the system when 
fli, and are restricted to a range appropriate to weak 
coupfing {ay,ar < 0.4). We label these the "dual reality" 
state and the "mixed reality" state. In the dual reality 
state, the oscillations of both pendula are uncorrelated 
and decrease in amplitude until both come to rest at the 
stable position = = 0. The two pendula behave as 
separate oscillators in the dual reality state, with reality 
and virtual reality interacting but coexisting individually. 
In the mixed reality state, the two pendula exhibit highly 
correlated stable, phase-locked periodic motion. In this 
mixed reality state of the system, the real pendulum and 
the virtual pendulum move together as one. These are 
illustrated in Fig. [21 

We model this coupled system by removing the real 
pendulum entirely and replacing it with a routine in the 
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code that separately integrates 

6sim + 2(39sim + ^siTU ~ Or fr COS 9 sim — 0, (3) 

an independent and equivalent equation of motion that 
represents the real pendulum. The (— ) sign with the 
Or term is necessary because the lever arm is mounted 
on a pivot that effectively reverses the direction of the 
movement of the actuator [see Fig.fljb)]. Since only the 
position of the real pendulum is measured in the experi- 
mental setup, the integration of Eq. ([3]) returns the posi- 
tion Osirn ■ The velocity is calculated using the position at 
the previous time step, and is evaluated analogously 
to Eq. ([2])- The virtual pendulum routine independently 
integrates 

4>sim + 2(3<j)sim + (wWo)^ Sin (j>sim + tty fy COS (j)sira = (4) 

at each time step. The only difference between the exper- 
imental system and the simulation system is that instead 
of measuring the position of the real pendulum, the sim- 
ulation integrates an equation of motion to calculate that 
position. The feedback works exactly the same way as 
before, except instead of an output voltage the program 
simply provides a feedback amplitude coefficient to the 
modeled real pendulum's integrator at each time step. 
This simulation can also run at real time, but suppress- 
ing the delay between integrations returns precisely the 
same results in a fraction of the time required. 

We now work with the reduced parameter space de- 
scribed by u) and a = y/ayUr, where a is the geomet- 
ric mean of the forcing amplitudes, a characterizes the 
strength of the bidirectional coupling; necessarily a — + 
as a„ — > or 0. For weak coupling, we have 

< a < 0.4. (D characterizes the quality of the model. 
For the virtual pendulum to be an accurate model of the 
real pendulum, Co has to be near 1. We find that there 
are two distinct regions in this parameter space, corre- 
sponding to two qualitatively different limiting behaviors 
of the system. These are depicted in Fig. [31 Region I cor- 
responds to the mixed reality state of the system. The 
oscillations are about the fixed point of each pendulum, 
and occur at frequencies close to the natural frequency 
of the real pendulum (see the top plot in Fig. [2]). Region 
II corresponds to the dual reality state of the system. In 
this region, both pendula initially oscillate but the sys- 
tem is unable to sustain this uncorrelated motion and 
loses kinetic energy until both pendula are at rest (see 
the bottom plot in Fig. [21). Region I has the Arnold 
tongue structure for mode-locked solutions in parameter 
space, as seen in Fig. [31 For each data set we wait until 

1 = 25 s, which is long enough for the transient dynamics 
to vanish, then we measure the maximum amplitude of 
the real pendulum X = max {9) over several oscillation 
cycles for each pair of parameters Q and a. We define 
X as the maximum amplitude of the real pendulum in 
the experimental system and Xgim as the maximum dis- 
placement of the simulated real pendulum. X and Xgi^ 
are plotted against (I) for a = 0.364 in Fig. [H This curve 



shows the phase transitions from Region II to Region I, 
then back to Region II. 
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FIG. 3: Limiting behavior phase transition diagram with 
Arnold tongue structure for the two parameters a and lu. The 
solid, dashed, and dotted lines indicate the critical points Uc 
in the experiment, simulation, and linear theory, respectively. 
These critical points are the boundary between Regions I and 
II in this parameter space. Region I corresponds to stable, 
phase-locked oscillations of the center of mass of each pendu- 
lum; this is the mixed reality state. Region II corresponds 
to both pendula ending at the stable equilibrium position 
(j) = = 0; this is the dual reality state. 

Dropping the sim subscripts, to linear order Eqs. ^ 
and (HI become 

9 + 2(39 + Lu^9 ~ ar4> = 0, (5) 
^ + 2(3<i) + 0j^ujl(j) + ay9 ^ 0. (6) 

By taking successive derivatives of Eq. ([5]) and substi- 
tuting these into the second derivative of Eq. this 
system can be written as the decoupled linear fourth- 
order differential equation (1 + a2)6'('') -I- 40^'^ + (4/3^ -|- 
-/^)9 + 2l3"f^9 + Q^u^9 = 0. 6i(4) = g^^) ^ ^ = at 
t — comprise the initial conditions, and we define 
7 = tt^o(l + w^). The general solution to this system is 
in the form 9{t) = '^i^'^'*, where the Ci are con- 

stants determined from the initial conditions and the 
Xi are the four solutions of the characteristic equation 
(1 + a2)r4 + 4:(3r^ + (4/32 + ^2)^2 _^ 2(3j^r + uo^uj^ = 0. 
If, for all i, Xi < 0, then we have 9{t) ^ as t ^ oo. 
However, even if for one eigenvalue we have Re(Ai) > 0, 
then the solution is no longer bounded: 9{t) ^ as 
t — > oo. For a given a, there are two values for Co such 
that Re(Ai) = for one eigenvalue. We define these to 
be the critical points Coc that mark the boundaries of the 
phase diagram in Fig. [3l As shown in this figure, the 
phase diagram of the linear system in Eqs. ^ and © 
closely matches that of the full system in Eqs. Q and (HI). 

There are differences between the experimental and 
simulation data. Region I in the simulation data extends 
further in the direction of small a than in the experimen- 
tal data (Fig. [HJ. Since the general shape of Region I is 
very similar for both sets of data, this onset of phase- 
locked oscillations for smaller a appears to reflect the 
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FIG. 4: This figure shows the response amphtudes X and 
Xsim versus a). The solid and dashed lines correspond to the 
experimental and simulation data, respectively. These are 
response curves for the system with d — 0.364. As indicated, 
for the experiment, the interval in uj where X > corresponds 
to Region 1, while the interval in ui where X = corresponds 
to Region II. Likewise, for the simulation, the interval in dj 
where Xsim > corresponds to Region 1, while the interval 
in ijj where Xsim ~ corresponds to Region II (see Fig. [S]). 

greater efficiency of the simulation. In the real system, 
there are inevitable small delays and noise in the elec- 
tronics and computer control, as well as additional fric- 
tion terms beyond the linear damping term in Eq. ([3]). 
As discussed above, more sophisticated models of fric- 
tion are necessary to reproduce pendulum dynamics ac- 
curately. Nonetheless, the virtual pendulum in Eq. ([T|) 
sufficiently models the real pendulum, giving rise to the 
mixed reality state when w « 1. 

This work presents experimental data from an interre- 
ality system comprising a virtual pendulum and its real- 
world counterpart. There is bidirectional, instantaneous 
coupling between the two pendula. We find that if the 
dynamics of the virtual system approximate those of the 
real one within a small tolerance, there is a phase transi- 
tion in the behavior of the system. The interreality sys- 
tem makes a transition from a dual reality state in which 
reality and virtual reality are uncorrelated to a mixed re- 
ality state in which reality and virtual reality are highly 
correlated. For both the experimental data and a numer- 



ical simulation, the region of mixed reality mode-locked 
solutions in parameter space is an Arnold tongue. While 
the appearance of an Arnold tongue is not surprising, 
it highlights two features of this interreality system that 
we would expect to be present in similar systems. The 
shape of the Arnold tongue indicates that with stronger 
coupling, the mixed reality states are accessible even with 
increased mismatch between the frequencies of the pen- 
dula. Also, the boundary of an Arnold tongue represent 
a discontinuous or sudden change change from a dual re- 
ality state to a mixed reality state, as seen in Fig. [H As 
virtual systems better approximate real ones, even weak 
coupling in other interreality systems may produce sud- 
den changes to mixed reality states. 

Forced systems tend to have the greatest response 
when the dynamics of the driving function match the 
dynamics of the system [l^l- This has been studied 
in de pth for various damped, coupled oscillator sys- 
tems One application of this phenomenon is res- 
onance spectroscopy, where a system with an unknown 
parameter b is driven with a forcing function that de- 
pends on this parameter. As b is varied, the system has 
the largest response when the dynamics of the forcing 
function match those of the system; thereby the value of 
b is identified p^ . To date, work in this area has predom- 
inantly focused on numerical analysis of coupled dilferen- 
tial equations that represent th e eq uations of motion for 
linear and nonlinear oscillators [ill. [l2l|. There has been 
some experimental work on non-sinusoidal driving of non- 
linear oscillators [lOj, but to the authors' knowledge there 
has been no work done on resonance spectroscopy via 
the instantaneous coupling of a virtual system to its real 
counterpart. It may be possible in this system or similar 
systems to use either the peak of the response curve or the 
boundary of the Arnold tongue to do system identifica- 
tion. Furthermore, it may be possible to experimentally 
observe universal routes to synchronization p^ . [l^ using 
interreality systems. We plan to examine these topics in 
the future. This work was supported by the National 
Science Foundation Grant No. NSF PHY 01-40179, NSF 
DMS 03-25939 ITR, and NSF DGE 03-38215. 
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